function model=update_alpha_2(model,n, MaxFun)


%options = optimset('LargeScale','on', 'MaxFunEvals', MaxFun, 'Display', 'off');

options = optimset('LargeScale','on','GradObj','on','Hessian','on', 'MaxFunEvals', MaxFun, 'Display', 'off');
x0      = model.alpha_2(n,:);
temp    = fminunc(@L_alpha_2, x0, options, model, n);

model.alpha_2(n,:) = temp;

end
